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Previous studies of phase mixing of ion cyclotron (IC), Alfvenic, waves in the collisionless regime 
have established the generation of parallel electric field and hence acceleration of electrons in the 
regions of transverse density inhomogeneity. However, outstanding issues were left open. Here 
we use 2.5D, relativistic, fully electromagnetic PIC (Particle-In-Cell) code and an analytic MHD 
(Magnetohydrodynamic) formulation, to establish the following points: (i) Using the generalised 
Ohm's law we find that the parallel electric field is supported mostly by the electron pressure 
tensor, with a smaller contribution from the electron inertia term, (ii) The generated parallel 
electric field and the fraction of accelerated electrons are independent of the IC wave frequency 
remaining at a level of six orders of magnitude larger than the Dreicer value and approximately 
20 % respectively. The generated parallel electric field and the fraction of accelerated electrons 
increase with the increase of IC wave amplitude. The generated parallel electric field seems to be 
independent of plasma beta, while the fraction of accelerated electrons strongly increases with the 
decrease of plasma beta (for plasma beta of 0.0001 the fraction of accelerated electrons can be as 
large as 47 %). (iii) In the collisionless regime IC wave dissipation length (that is defined as the 
distance over which the wave damps) variation with the driving frequency shows a deviation from 
the analytical MHD result, which we attribute to a possible frequency dependence of the effective 
resistivity, (iv) Effective anomalous resistivity, inferred from our numerical simulations, is at least 
four orders of magnitude larger than the classical Spitzer value. 

PACS numbers: 52.35.Hr; 96.50.Tf; 04.30.Nk; 96.50.Ci; 96.60.P-;52.65.Rr; 96.60.pf; 96.60.qe 



I. INTRODUCTION 

Phase mixing is a mechanism of enhanced dissipation 
of Alfven waves due to inhomogeneity of Alfven speed in a 
direction transverse to a local magnetic field. This mech- 
anism originally was studied in the fusion and laboratory 
plasma context by a number of authors P, [2j_^, U, d, Q 
and subsequently applied to the solar corona 3 ■ Most of 
the large amount of work done in the field of phase mix- 
ing was in the resistive MHD (Magnetohydrodynamic) 
regime. Recently, a few works looked at the same mech- 
anism in the collisionless regime in the context of Earth 
magnetosphere [Qj and solar corona [l^, El ■ The main 
findings of these works include the generation of electric 
field that is parallel to the ambient magnetic field in the 
regions of transverse density inhomogeneity, as well as 
associated electron acceleration. It should be mentioned 
that these studies considered circularly polarised ion cy- 
clotron (IC) waves which in the low frequency regime 
become Alfven waves. We use terms Alfven or IC inter- 
changeably, but reader should bear in mind we always 
refer to waves with frequencies < ojd (with ojd being ion 
cyclotron frequency). The exact mechanism of genera- 
tion of the parallel electric field has stimulated a debate 
[13, , ^'^'i even MHD regime option was explored [ij] . 
Continuing this investigation, here we apply technique 
used in the collisionless reconnection [is., i6]. Namely, 
in section 3.1 we use generalised Ohm's law to find out 
which term generates the parallel electric field. 

Solar flare observations [l3l trigger one's interest in 
how effectively plasma particles are accelerated. Hence, 
in section 3.2 we look into how the generated parallel 



electric field and the fraction of accelerated electrons de- 
pend on model parameters such as IC wave frequency, 
amplitude, and plasma beta. 

Ref. [71] provides a simple analytical expression how 
Alfven wave amplitude should decay in space due to the 
phase mixing. Despite the fact that their formula is de- 
rived in the resistive MHD regime, we still apply it to 
our collisionless, kinetic simulation and see what does 
the comparison yield (Section 3.3). This is done in the 
light of previous results of Tsiklauri et al. who es- 
tablished that in the collisionless, kinetic regime Alfven 
wave amplitude in the density gradient regions decays 
with distance (from where it is driven) according to col- 
lisional MHD formula of Heyvaerts and Priest [7| . Here 
we stretch the MHD-kinetic analogy further to test 
dependence under the exponent. 

In Sect 3.4 we estimate the effective "resistivity" (again 
the spirit of MHD-kinetic analogy) . The quotation marks 
are needed to signify that PIC (Particle-In-Cell) simula- 
tion code is collisionless and hence no resistive effects ex- 
ist as such. However, scattering of particles by magnetic 
fields plays effective role of collisions. 



II. SIMULATION MODEL 

In our numerical simulations we use 2.5D, relativis- 
tic, fully electromagnetic PIC code. The initial condi- 
tions, basic parameters and boundary conditions arc ex- 
actly the same as in the previous work by Tsiklauri ct al. 
p^ . In particular, the uniform magnetic field is in x- 
direction, the transverse density inhomogeneity is across 
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y-direction. z is the spatially ignorable coordinate, i.e. 
djdz = 0. However, we retain all three components of 
velocity, v^-, Vy and Vz- The system size without ghost 
cells in two dimensions is Lx — 5000A and Ly = 200A 
where A(= 1) is the simulation grid size, corresponding 
to the electron Debye length, A_d = vte/i^pe — lA (vte is 
electron thermal velocity and ujpe is electron plasma fre- 
quency). The total particle number is 2.39 x 10* electron- 
ion pairs. The ion to electron mass ratio is mi/nie — 16 
due to the limitation of computer resources and speed. 
Tsiklauri et al. [TO] use the fixed driving wave amplitude 
SB/Bq = 0.05, and plasma j3 — 0.02. The dimensionless 
ion and electron density inhomogeneity is described by 
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These are normalised to some background constant value 
(no). Here all plasma parameters are quoted at the 
boundary, away from the density inhomogeneity region. 
In the central region (across y-coordinate) , the density 
is smoothly enhanced by a factor of 4, and there are the 
strongest density gradients having a width of about 51A 
around the locations y = 51. 5A and y — 148. 5A. Below, 
in Eqs.(4) and (5) we shall be using y ~ 51.5 — 51/2 ~ 
26A and y = 51.5 -f 51/2 = 77 A as the boundaries of 
one of the transverse density gradients. The background 
temperature of ions and electrons, and their thermal ve- 
locities are varied accordingly 

T,{y)/To^Te{y)/To = F{y)-\ 

Vth,i/ViO = Vth.e/VeO = F(?/)~^^^, 

such that the thermal pressure remains constant. Since 
the background magnetic field along the a;-coordinate is 
also constant, the total pressure remains constant too. 
Then we impose a current of the following form 

dtEy = -Jo siniudt) (1 - exp [-(t/to)'^]) , 

dtEz = -Jocos{ujdt) (1 - exp [-{t/to)^]) . 

In Ref. [l3] the driving frequency was fixed at ivd — 

0. 3(jJci- Here, we also use driving frequencies that satisfy 
(jjd < ujci so that no significant ion-cyclotron resonant 
damping takes place, dt denotes the time derivative, 

is the onset time of the driver, which was fixed at SO/wpe 

1. e. 3.125/cL;ci. This means that the driver onset time is 
about 3 ion-cyclotron periods. Imposing such a current 
on the system results in the generation of left circularly 
polarised IC (Alfvenic) wave, which is driven at the left 
boundary of simulation box and has spatial driver width 
of lA. The wave propagates along a:-coordinate and gen- 
erates the parallel electric field in the density gradient 
regions (for more details see Tsiklauri et al. [ifll)- The 
parameters used are commensurate to what is seen in so- 
lar corona by e.g. Doppler broadening of emission lines. 



The observed values of the Alfven waves at heights of 
R = 1.04i?sun = 28 Mm are about 50 km s~^ (see e.g. 
[Ts!]). which for a typical Alfven speed of 1000 km s~^ 
makes SB/Bq equal to 0.05. 

As one of our goals is to investigate the parame- 
ter space of the problem, the following range of phys- 
ical parameters was used: we vary the frequency, am- 
plitude of IC wave or the plasma beta. The driving 
wave frequency LOd/uJci ranged from 0.1 to 0.5 with a 
step of 0.1. The wave amplitude SB/Bq was also set 
at 0.01,0.05,0.10,0.15,0.20 and 0.25. The plasma beta 
varies from f3 — 10^^ to 10^^ by controlling the elec- 
tron thermal velocity. Specifically, the plasma beta 
was set at 0.0001,0.0003,0.0010,0.0030,0.0100,0.0200 
and 0.0300. Each beta value corresponds to Vte/c = 
0.007,0.012,0.022,0.039,0.071,0.100 and 0.122, respec- 
tively, where vte is taken from the low density region and 
c is speed of light. Each numerical run (each data point 
in subsequent Figs. (2)-(5) typically takes about 7-10 
days on 64 parallel processors. 



III. SIMULATION RESULTS 

A. Source of the parallel electric field 

In order to understand details of the parallel electric 
field generation process, we now focus on the question: 
which term in the generalised Ohm's law is responsible 
for the generation of the parallel electric field? The gen- 
eralised Ohm's law can be written as 



E = -v,xB £_^^+^/^.VK, 1 

UpC e at 



where E and B are electric and magnetic fields, v is 
plasma velocity, P is pressure tensor (3x3 matrix), n 
is plasma number density, m is mass and e is electric 
charge. The subscript e refers to an electron. Normal- 
ising space coordinate by = 5000 A, fluid velocity by 
Alfven speed (at the lowest density region) vao: time by 
Alfven transit time ta{= L/vao), magnetic field by Bq^ 
electric field by vaoBq^ number density by uq and pres- 
sure tensor by Bq/moj a dimensionless version of Eq. ^ 
can be obtained 



-> V ■ P Tfi Qv 

E=-VpXB- d, - d, — {—^ + {Ve ■ W)Ve) , (2) 

rip m, ot 



where di is the normalised ion skin depth {di = c/ujjnL). 
Note that strictly speaking we should have used tildes in 
Eq. ([2]) to denote dimensionless quantities, but we omit 
them for brevity. 

Let us focus on the parallel electric field, E^, which is 
generated in the density gradient regions by phase mix- 
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FIG. 1: Different term profiles along i-coordinate 
from the generalised Ohm's law: (a) Ex (normalised to 
VAoBo), (b) -{di/n)dPexx/dx, (c) -{di/n)dPeyx/dy, and (d) 
—di{me/mi)dvex/dt in the maximum density gradient re- 
gion y — 0.0103 at t — 0.4375 for the case, uud/ud = 0.3, 
SB /Bo — 0.25 and (3 — 0.02. The time derivative term in 
(d) is calculated from t = 0.4370 to 0.4375 (time interval cor- 
responds to the inverse of electron plasma frequency 
Refer to text for the normalisation used. 



ing. It is given by, 



Ex = -{VeyB^ - VezBy) - di- 



1 /BP BP 

^ I exx , eyx 



n \ Bx dy 



me / BVex , dVex , dVex , 
-di [^:r—+Vex— ^Vey— ), (3) 



rrii \ dt 



Bx 



dy 



where 8/Bz = is assumed because of spatially 2 dimen- 
sional model. 

Fig. [1] shows the different term profiles (along the uni- 



form magnetic field in cc-direction) in the generalised 
Ohm's law (Eq. §1): (a) Ex, (b) -{dJn)8Pexx/dx, 
(c) -{di/n)BPeyx/dy, and (d) ~di{me/mi)8vex/dt in 
the maximum density gradient region y — 0.0103 at 
t = 0.4375 for the case, = 0.3, SB/Bq = 0.25 

and P = 0.02. The time derivative term in Fig. [Hd) is 
calculated from t = 0.4370 to 0.4375 (time interval corre- 
sponds to the inverse of electron plasma frequency ^^pe)- 
The other terms in right-hand side of Eq. ([3]) are neg- 
ligible. In Fig. [TJa), the parallel electric field is clearly 
observed in the density gradient regions where phase mix- 
ing can occur. It should be noted that no parallel electric 
field is seen away from the density gradient regions. Ac- 
cording to the generalised Ohm's law, there has to be a 
term on the right-hand side of Eq. ^ that supports this 
electric field. 

By comparing Fig. [Ija) to Fig. [Hb-d) it is clear that 
the parallel electric field (Fig. [Ija)) is supported mostly 
by the electron pressure tensor (Fig.[T](b)), with a smaller 
contribution from the electron inertia term (Fig. [IKd)). 
It is interesting to note that previous results on collision- 
less reconnection both in tearing unstable Harris current 
sheet [H [H, [131 and stressed X-point collapse [H, lU 
have shown that the term in the generalised Ohm's law 
that is responsible for breaking the frozen-in condition, 
i.e. that supports out-of-plane electric field is electron 
pressure tensor. Here the electron pressure tensor sup- 
ports (generates) the parallel electric field. Thus, one 
should note a universal importance of the electron pres- 
sure tensor in relation to supporting the electric fields in 
coUisionless plasmas. 



B. Parametric study of the generated parallel 
electric field and the fraction of accelerated electrons 

To estimate how efficiently the parallel electric field is 
generated by phase mixing, as a function of model pa- 
rameters, we introduce the average of the absolute value 
of the parallel electric field in the density gradient region: 

r-L^ p77A 

(4) 



E* 



1 



LxLy 



\Exix,y)\ 



x=0 Jy=26A 



En 



dxdy, 



where Eq ~ mi^cLOp,,/ g. Note that in what follows the 
normalisation of the electric field is different from Sec- 
tion 3.1 where, usual for the generalised Ohm's law, 
"fluid" normalisation (vaoBq) is used. Normalisation 
Eq — meOjUpe/e is usually referred to as "electrostatic". 
By using the definition given by Eq.jjl]) we can evaluate 
quantitatively the electric field generated in the density 
gradient region. Although there are two density regions 
in our simulation box because of the periodic boundary 
condition, we focus on the lower one (the physics of the 
upper and the lower regions is essentially the same). The 
range from y — 26A to 77A indicates the density gradi- 
ent. See Tsiklauri et al. [lH for details. 

Also, in order to investigate the fraction of acceler- 
ated electrons by the generated parallel electric field. 
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we use particle data in the density gradient at lower 
side (26 A < y < 77A). We count the number of elec- 
trons with velocities greater than the thermal velocity 
{vte < Vx < c) in the electron velocity distribution 
function, in the x-direction, at the final time snapshot 
ujcit — 54.69, and divide this by the total number of par- 
ticles (with < Vx < c) in the same domain: 



(q) 



N 



r 

JVx=V 



rL^ r77A 

jx=0 Jy=26A 



f{vx)dvxdxdy 



r r^" r^^"^ 

Jv^=0 Jx=0 Jy=2eA 



f{vx)dvxdxdy 



(5) 



Here it was to enough to integrate only positive region in 
this distribution because electron acceleration was sym- 
metrical in the x-direction. Note that the initial veloc- 
ity distribution function is nearly Maxwellian. To main- 
tain the balance of the total kinetic pressure throughout 
the system, the particle thermal velocity in the dense 
plasma region {y — lOOA) is lower than the outside re- 
gion {y = OA or 200 A) (see for details Fig. (4) from Tsik- 
lauri et al. However, in order to use Eq. we 

had to estimate an appropriate velocity corresponding to 
the thermal velocity in the Maxwellian. Fortunately, ini- 
tial electron velocity distribution function did not deviate 
much from the exact Maxwellian. Therefore, we adopted 
the standard thermal velocity which was set to 36.8% of 
f{v — 0). Recall that f{v — vte) = noexp{—v^/vf^) — 
noexp(— 1.0) = 0.368no, where uq is the peak number at 
V = in Maxwellian. 

The choice of diagnostic for characterizing the de- 
gree of particle acceleration needs a further clarification. 
Eq.(5) provides number of electrons with speed Vx ex- 
ceeding the thermal speed Vte as a fraction of total dis- 
tribution value. However, in a Maxwellian plasma 16% 
of the electrons have v > Vte (for vte = y^feBT^/mT). 
Thus, it may be tempting to either: (i) subtract this 
16% from our diagnostic in Eq.(5) or (ii) try to fit a 
Maxwellian to the electron distribution at the final sim- 
ulation time step (after all the acceleration has taken 
place) and then count the number of electrons that have 
speeds Vx > vte- Our motivation to keep the definition 
given by Eq.(5) is two- fold: (i) In context of the solar flare 
observations in X-rays, one always infers the integrated 
spectrum averaged over some volume i.e. < f{v)nV > 
(where n is electron number density) and it is impossible 
to subtract the above mentioned 16% without introduc- 
ing additional uncertainties due to the unknown density 
and poorly known emitting volume, due to a line of sight 
effect (Dr. E. Kontar of Universityof Glasgow, private 
communication), see also Refs. p2l.l23l|. Hence the defini- 
tion given by Eq.(5) is more appropriate for comparison 
of theory with the observations. Also, this is particularly 
timely because the acceleration of electrons by Alfvenic 
waves in flares have been recently studied [2J]. (ii) Al- 
though the electron distribution function at the initial 
time step is nearly Maxwellian (despite density and tem- 
perature transverse inhomogeneities), at the final stages 
of the simulation the deviations from the Maxwellian 
form are large and hence a fit to a Maxwellian, in order 
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FIG. 2: (a) E* /Eq vs. ujd/i^ci and (b) A''/A''o vs. LUd/uid for 
SB /Bo = 0.05 and /? = 0.02. 



to calculate Vte, and in turn to count the super-thermal 
particles above that value {vx > vte), is impractical. 

Fig. ([2) shows the generated parallel electric field, 
E* /Eq, and the fraction of accelerated electrons, N/Nq 
as function of driving frequency of the IC wave (nor- 
malised to LOci). The wave frequency ujd/ujci was set at 
0.1, 0.2, 0.3, 0.4 and 0.5. The wave amplitude and plasma 
beta were fixed at 5B/Bq = 0.05 and (3 = 0.02. We 
gather from Fig. ^ that the generated parallel electric 
field and the fraction of accelerated electrons are indepen- 
dent of the IC wave frequency remaining at a level of six 
orders of magnitude larger than the Dreicer value and ap- 
proximately 20 % respectively. This conclusion is based 
on the following estimate E*/Eo « 0.0018, i.e. E* = 
5473.4337 V m^^ (note that here solar coronal number 
density of n = 10^^ m~'^ was used). The Dreicer electric 
field (which is associated with the particle acceleration 
run-away regime [111), Ed — (ne^ In A)/(47regfc_BT), for 
T = 1 MK (and hence InA = 18.095) is Ed = 0.00547 V 
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J: (a) E'/Eo vs. SB /Bo and (b) N/Nq vs. SB/Bq for 

c^dMi = 0.3 and /3 = 0.02. 



FIG. 3: 



FIG. 4: (a) E*/Eo vs. /3 and (b) iV/7Vo vs. (3 for Wd/a;^ = 0.3 
and 5B/Bo = 0.05. 



m-^ Thus, E*/Ed = 10^. This result should be taken 
with caution because it is obtained with the ion-electron 
mass ratio of 16. As shown by Tsiklauri [l^ (see their 
Figure 7) that attained amplitude of the generated par- 
allel electric field scales strongly as oc 1/ {rni/m^). 

Next, we explore how the generated parallel electric 
field and the fraction of accelerated electrons depend on 
driving IC wave amplitude. The latter, SB/Bq, was set at 
0.01,0.05,0.10,0.15,0.20 and 0.25. The wave frequency 
and plasma beta were fixed at Lo^/uj^i = 0.3 and ji = 0.02. 
We gather from Fig.[3]that the generated parallel electric 
field and the fraction of accelerated electrons increase 
with the increase of IC wave amplitude. This seems as 
a reasonable result because larger amplitude waves have 
more energy to give to electrons. Also, non-linear effects 
would be progressively important. 

We also explored the plasma beta de- 
pendence. The plasma beta was set at 
0.0001,0.0003,0.0010,0.0030,0.0100,0.0200 and 0.0300. 



The wave frequency and amplitude were fixed at 
^dl^ci = 0.3 and 5B/Bq — 0.05. The plasma beta 
is defined as = 2fiop/B'^ — 2{vte/c)'^ /{uJce/t^pe)^ ■ 
We altered the plasma beta by changing the electron 
thermal velocity, affecting the plasma kinetic pressure. 
Therefore, in this simulation, magnetic pressure is kept 
constant while plasma kinetic pressure varies. Fig. [H^a) 
does not show any correlation between plasma beta 
and the parallel electric field generation. Incidentally, 
Tsiklauri et al. [26] investigated plasma beta dependence 
of the fast magnetosonic wave amplitude, which is gen- 
erated in a transversely inhomogeneous medium when 
an Alfvenic pulse is launched (using MHD numerical 
simulation). According to Fig. 9(b) in Tsiklauri et al. 
[26| . the maximum fast magnetosonic wave amplitude 
also does not depend on plasma beta. Tsiklauri [l3| 
alluded to the relation between the non-linear fast 
magnetosonic wave and parallel electric field generation. 
Hence it is not surprising that in our Fig. IDJa) we do 
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FIG. 5: A log-log plot of the dissipation length Ld/{c/upi) 
vs. driving IC wave frequency LOd/iOd (open symbols). The 
solid line is the least squares fit with a slope of —1.10. 



not see plasma beta dependence. What is surprising, 
at the first glance, in Fig. [H^b), is that the fraction 
of accelerated particles strongly depends on beta. In 
particular, a decrease in beta (for (3 — 0.0001) yields 
as large percentage as N/Nq — 0.472 ^ 47%! One 
can conjecture that this is due to the fact the in the 
case of small plasma beta, magnetic effects dominate 
over thermal ones, and because IC wave is essentially a 
magnetic-type perturbation, electrons respond better to 
the wave influence and accelerate more efficiently. 



In particular, one can e.g. check the scaling in the 
Eq.®. For this purpose, after simple algebra applied to 
Bz oc exp[— (x/Ld)'^], for fixed rj^dyVA and va, one can 
obtain the following scaling logj^Q Ld — —2/3 log^Q Ud + C, 
where C is some constant. Hence the slope of a log-log 
graph of dissipation length Ld versus ujd is expected to 
be —2/3 = —0.67. Fig. [5] presents the dissipation length 
(distance over which the wave damps) dependence on 
the IC wave frequency. The values in Fig. [5] were ob- 
tained from line data of magnetic field Bz in the maxi- 
mum plasma gradient Y/(c/u>pe) = 14.8 at LUdt — 82.0 for 
^dl^ci = 0.2, cjcji = 54.7 for lOdl^d = 0.3, lOdi = 41.0 
for LOd/iOci = 0.4 and ojcit — 32.8 for ojd/iOci = 0.5, re- 
spectively. The reason for the different snapshot times is 
that we only consider well-developed phase mixing, i.e. 
when IC wave is fully damped. As in left panel of Fig. 6 
in Tsiklauri et al. [l^) we fit Bz{x) to oc exp {~{x/ LdY) 
and obtain empirical dissipation length Ld- we gather 
from Fig.[5]that the slope is —1.10 contrary to the above 
prediction of —2/3 = —0.67. 

To address the inconsistency we conjecture that the 
resistivity might be variable. One can estimate the re- 
sistivity for each case of driving wave frequency consid- 
ered, by calculating rj = (!>v\/ {uj'^{dyVA)'^ L\). We reit- 
erate that strictly speaking PIC simulation code is col- 
lisionless and hence no resistive effects exist. However, 
scattering of particles by magnetic fields plays effective 
role of collisions. When resistivity is mentioned we refer 
to "effective" resistivity. Normalising the frequency by 
IC frequency uJd, length by ion skin depth c/upi and ve- 
locity by speed of light c, the dimensionless resistivity is 
given by. 



C. Amplitude decay law in the kinetic regime 

Tsiklauri et al. [13] established that in the collision- 
less, kinetic regime Alfven (IC) wave amplitude in den- 
sity gradient region decays with distance (from where it 
is driven) according to collisional MHD formula of Hey- 
vaerts and Priest [3| 



Bz exp 



iqUjl{dyVAf 

6v\ 



(6) 



where rj is resistivity (divided by /ig, i.e. by 77 we mean 
l/(CT/io), is (driving) wave frequency, va is Alfven ve- 
locity and X is the axis which AW propagates. Recall 
that plasma density inhomogeneity is a function of y and 
wave propagates in the x direction. In particular, (see for 
details left panel of Fig. 6 in Tsiklauri et al. |10|]) it was 
shown that, at a fixed time instance corresponding to 
well-developed phase mixing, Bz{x) oc exp (— (x/L^)'^) . 
In other words, one can define an empirical dissipa- 
tion length, Ld, according to collisional MHD formula 
(Eq.®) in the coUisionless, kinetic regime. Applicability 
of the MHD formula in the kinetic regime was rather sur- 
prising. We now turn back to the previous results in order 
to see how far this MHD-kinetic analogy can be stretched. 



V 



27 V w, 



\ 



1 + 3 exp 



y-yc 

L 



exp 



y-yc 



y-yc 

L 



(7) 



We can now put in the known parameters L/(c/LUpi) — 
1.25, va/c = 0.25 and {y - yc)/L = (148.5 - 100)/50 = 
0.97 into Eq. Here L and j/c are the width and the 
centre of plasma density gradient, respectively. Therefore 
the dimensionless resistivity can be estimated using 



V 



= 1.97 X 10"^ — 



^d 



c/ LUp, 



(8) 



As the driving wave frequency and the dissipation length 
for each case are given by choice and empirically, re- 
spectively, one can obtain the resistivity by substituting 
above values into Eq. ([8]). The result is r]/{c'^ujpi/ujd) = 
5.88 X lO"'^ (for Wd/wc» = 0.2), 9.21 x 10"^ (for uJd/(^cz = 
0.3), 1.28 x 10-5 (for Wrf/w^ = 0.4) and 2.02 x lO''^ 
(for LUd/ujci — 0.5). Thus our initial conjecture that the 
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effective resistivity depends on the driving IC wave fre- 
quency turns out to be correct. But the main conclusion 
of this analysis is that despite collisional MHD scaling 
Eq.® being applicable to the coUisionless, kinetic regime 
of phase mixing, i.e. Bz{x) oc exp { — {x/ Ld)^) scaling 
holds, stretching the MHD-kinetic analogy further to 
dependence under the exponent is not valid (due to the 
effective resistivity being a function of w^). 

D. Effective anomalous resistivity 

Issue of anomalous resistivity is central for many space 
and laboratory plasma applications. It can facilitate fast 
magnetic reconnection via Petschek type mechanism (if rj 
is not spatially uniform), or have significant implications 
for wave heating models of solar corona where normal 
Spitzer resistivity is too small to produce any sizable ef- 
fect. Ref.[27| presented plasma resistivity measurements 
in the reconnection current sheet of Magnetic Reconnec- 
tion Experiment (MRX) (see ref. for details of the 
experimental setup). They established that in the coUi- 
sionless regime measured resistivity values can be more 
than an order of magnitude larger than the Spitzer value 

Let us apply our PIC simulation results to see if there is 
any evidence for the anomalous effective resistivity. We 
fix physical parameters corresponding to solar coronal 
plasmas: B — O.OIT and plasma number density hq — 
2 X lO^^TO^^^ i.e. LOci = eB/{l6me) = 1.10 x lO^rad/s 
and LUp, = ^noe2/(16meeo) = 6.30 x lO^rad/s. Eq. (j8]) 
can be rewritten as: 

rj = 9.2AxW(^y' (4^y\ (9) 

\UJci J \c/UJpiJ 

Similarly to the previous calculation in Eq. ([8]), now us- 
ing Eq. jH) we obtain -q = 2.76 x lO" (for cod/c^d = 0.2), 
4.32x10"* (forwd/wci = 0.3), 6.02x10'' (forwd^i = 0.4) 
and 9.46 x lO'' (for uJd/oJd = 0.5). Here units of the resis- 
tivity are m^ sec"-'^. We gather that all values are in the 
range of « lO'' — 10^ m^ scc^^. Spitzer resistivity (nor- 
malised to /io) for the above parameters and T = 1 MK 
is 1.83 m^ sec^*. Thus, we conclude that our numerical 
simulations provide effective resistivity values of w lO'' 
times larger than Spitzer value, which is indicative of the 
anomalous resistivity. It should be mentioned that these 
results were obtained for the ion-to-electron mass ratio 
of 16. Clearly one would expect some dependence of the 
obtained effective resistivity on the mass ratio. Thus, the 
obtained results should be taken with caution. 



IV. CONCLUSIONS 

Let us summarise the above findings: 



We used the generalised Ohm's law and found that the 
parallel electric field, which is generated by propagation 
of IC (Alfvenic) wave in a transversely inhomogeneous 
plasma, is supported mostly by the electron pressure ten- 
sor, with a smaller contribution from the electron iner- 
tia term. Surprisingly, this result resembles closely to 
the previous results on coUisionless reconnection both in 
tearing unstable Harris current sheet [l^, [H, H^l and 
stressed X-point collapse [H, [2lj. However, in the lat- 
ter two cases, the generated electric field is in the plane 
perpendicular to the magnetic field. Thus, a universal 
importance of the electron pressure tensor in relation 
to supporting the electric fields in coUisionless plasmas 
should be noted. 

We explored physical parameter space of the problem 
with regards to the efficiency of generation of parallel 
electric field and acceleration of electrons. We found that 
the generated parallel electric field and the fraction of 
accelerated electrons are independent of the IC wave fre- 
quency staying at a level that is 10^ times larger than the 
Dreicer value and approximately 20 % respectively. The 
generated parallel electric field and the fraction of ac- 
celerated electrons increase with the increase of IC wave 
amplitude. The generated parallel electric field seems 
to be independent of plasma-/?. However, the fraction of 
accelerated electrons strongly increases with the decrease 
of plasma-/?, e.g. for plasma /3 = 0.0001 the fraction of 
accelerated electrons can be as large as 47 %. 

Previously it was established that in the coUision- 
less, kinetic regime phase-mixed Alfven (IC) wave am- 
plitude damps with distance of propagation according 
to cx exp[—{x/Ld)^] which resembles closelyto col- 
lisional MHD result of Heyvaerts and Priest [3]. We 
tried to stretch this analogy further by investigating how 
the dissipation length Ld scales with the IC driving fre- 
quency. We found that the scaling is different from the 
MHD result. We have shown that this discrepancy can 
be attributed to the frequency dependence of the effective 
resistivity. 

We also found that the effective resistivity, albeit for 
unrealistic mass ratio, still is as large as 10"* times the 
classical Spitzer value. 
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